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TECHNICAL MEMORANDUM 


CALIBRATION OF THE ISEE PLASMA COMPOSITION EXPERIMENT 
I. APERTURE POINTING AND APERTURE ACCEPTANCE ANGLE 


Figure 1 is a drawing of the ISEE spacecraft. The spacecraft is a 16-sided 
polygon with booms. The +Y-axis is along the ACS extention and emerges from the 
spacecraft body between Faces 1 and 16. The +X-axis emerges from the spacecraft 
body between Faces 12 and 13. The Plasma Composition Experiment is actually two 
independent and almost identical instruments; Head A and Head B. Both heads are 
in the same face mounted one above the other. Head A is on top and Head B is 
mounted below but is upside down with respect to Head A. Although each head looks 
radially outward, they are canted slightly such that they do not look perpendicular 
to the spacecraft spin axis (+Z). The center line of Head A is directed 95 deg from 
the spin axis, the center line of Head B is directed 85 deg from the spin axis. 

The ISEE orbit is inclined 30 deg to the Earth's equator. The spacecraft's spin 
axis is actively controlled to within about a degree of perpendicular to the ecliptic 
plane. Hence, the two heads are looking approximately in the ecliptic plane, but the 
relation of the pointing of the heads with respect to the velocity vector is more 
complicated. Figure 2 is a predicted plot of the angle between the perpendicular to 
the spin axis and the velocity vector as a function of distance from the center of the 
Earth during the first four years of the spacecraft's lifetime. 



Figure 1. ISEE-1 satellite with principal appendages. 



VELOCITY VECTOR ECLIPTIC LATITUDE PROFILE 
SHOWN AT 600 DAY INTERVALS 
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Figure 2. Synopsis of the ISEE-1 satellite orbit. 


Figure 3 is a diagram of the coordinate system used to describe the geometry 
of the experiment. The various angles are defined as follows: 

0 = THETAS = the angle between the Z-axis and the velocity axis 

o 

0p = THETAD = the angle between the Z-axis and the center of the detector 
aperture 

= PHID = the angle between the X-axis and the center of the detector 
aperture 

NOTE: The coordinate system is defined such that the spacecraft velocity 

lies in the X-Z plane. 

a = ALPHA = the azimuthial acceptance angle 
3 = BETA = the polar acceptance angle 


Each head of the Plasma Composition Experiment consists of a rectangular aper- 
ture (9.4 x 2.2 cm) which is mounted flush with the spacecraft surface. Following 
the aperture is a collimating section which, ideally, establishes the angular response 
characteristics of the instrument purely by its geometry. The aperture acceptance 
angle is nominally a = 12 deg, 3 = 45 deg. 

Following this section is a region of electrostatic acceleration in which the 
incoming ions are accelerated to approximately 3000 electron volts. Subsequently, this 
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Figure 3. Definition of satellite coordinate system. 

ion beam is energy analyzed and split. A portion is detected by channel electron 
multipliers and a portion is input to the mass spectrometer, mass analyzed, and then 
detected by a Johnston multiplier. 

Two independent attempts have been made to analytically model the instrument. 
The first is a simple first order analysis restricted to an analysis of the collimating 
section of the instrument using a Maxwellian input. The second method was an 
elegant modeling of the trajectories of individual particles through the fields of the 
electrostatic analyzer section of the instrument. The first order approach is pre- 
sented first, followed by the Head A calibration results, and then the ray tracing 
effort. 

To a first approximation, it can be assumed that once the incoming ions are 
accelerated by the 3000 V potential, they "fly straight" through the subsequent 
analyzers and detectors. This assumption has the effect of "imaging" the detector 
after the energy analyzer and the detector after the mass analyzer forward to the 
plane at which the ions are accelerated. The resulting configuration is shown in 
Figure 4(a,b,c). (See also Figure 22 in the Ray Tracing section.) Figure 4a is a 
vertical slice through the instrument in the plane of the Z-axis and the detector 
normal. The aperture is to the left (0.0944 m) ; the detectors are arrayed to the 
right and labeled MD for the mass analyzer detector, ED 2 for the high gain energy 
analyzer detector, and EDI for the low gain energy analyzer detector. It is to be 
noted that the switch between high and low gain in the energy analyzer detector is 
accomplished by actually physically switching between a pair of detectors with differ- 
ent effective areas; whereas the change in gain at the mass analyzer detector is 
accomplished by changing the bias on the plates of a single detector. 
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Figure 4a. PCE diagram, side view, with emphasis on electrostatic analyzer. 
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Figure 4b. PCE diagram, side view and top view, with emphasis 

on optical path dimensions. 
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Figure 4e. PCE diagram, energy detector apertures. 

Figure 4: The instrument consists of: 

1) Collimater 

2) RPA: 

a) Ground grid 

b) Retarding grid 

c) Ground grid 

3) Pre- acceleration region — 3 kV (actually varied in 25 V steps around 
-2950 V — see Energy Passband Section) 

4) ESA 

5) Magnets 


5 




Note the relative placement of the three detectors. MD is in the most advanta- 
geous position, directly centered behind the limiting aperture opening. The angle 
subtended between the center of this detector and the aperture (angle 6) is approxi- 
mately 95 deg with the center line of this angle down the "bore- sight" of the instru- 
ment. EDI is offset from the center line. It can be calculated that the total angle 
(6) subtended between EDI and the aperture is still approximately 45 deg but from 
the geometrical considerations, the center line of the angle has been shifted approxi- 
mately 12 deg off the center line of the instrument. 

ED 2 is in the most compromised location. Its surface is actually partially 
masked from the aperture by structure internal to the instrument. The total angle 
subtended by the aperture at the center of this detector remains close to 45 deg 
(6 = 43 deg is actually calculated) , but the center line of this detector's viewing is 
shifted 16.5 deg above the instrument's center line. 

Two further complications should be considered. The first is the angle of 
incidence of the particles when they strike the detector. Particles which enter the 
aperture at the upper edge and strike ED 2 and are counted are arriving at the detec- 
tor at a 38 deg angle of incidence. Particles which enter the lower edge of the 
aperture strike the detector at near normal incidence. The second effect is due to 
the finite size of the detectors. Each point on the detector face actually has a 
slightly different viewing angle through the aperture. This effect would be most 
pronounced for MD since it is physically the largest, although it could be argued that 
it might also influence the energy detectors. 

To check these two effects, a calculation was made in which the face of each 
detector was subdivided into a square array of elements. The angle subtended 
between each element and the aperture was calculated and the response of this 
"elemental," limited aperture to a drifting Maxwellian was calculated. The results 
from each of the elements were then summed. The calculation was repeated for 
different THETAD angles with PHID = 0 to form a response versus angle curve. 

Of the two effects mentioned, only the cosine effect was significant, and this only 
for ED 2. All detectors were physically small enough to be considered, for all 
practical purposes, point detectors. The results of the calculation for ED 2 are 
shown in Figure 5. Three of the curves in this plot are calculated 6 response 
curves, the fourth curve is calibration data taken utilizing a 19 eV beam for Head B. 

Two of the calculations were made at high Mach numbers to simulate the calibra- 
tion beam. The third curve was calculated at a low Mach number to demonstrate that 
the effects which are being discussed, and even the relative importance of the 
absolute pointing, tend to lose significance at the lower Mach numbers. 

It is interesting to note that although the center of the full- width-half-max of 
the curve is indeed centered approximately 16 deg off the center line of the instru- 
ment (90 deg point on the plot), the ED 2 response function tends to be higher near 
the center line than away from it. To first order, this phenomenon would tend to 
shift the ED 2 average response toward the instrument's center line by a few degrees, 
bringing it more into line with EDI, where the effect of an asymmetric peaking was 
less pronounced. 

All of the above is valid only insofar as it is supported by calibration data; 
which in the case ED2, Head B, it was not. The calibration curve on the plot shows 
that the measured response of the instrument appeared to be down substantially near 
the instrument center line. This points to the possibility that the detector is actually 
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Figure 5. Detector angular response calibration, Head B, 
energy detector 2, and model. 

cut off from the aperture by the internal structure by more than was anticipated or 
that strong fields inside the collimating section are playing an important role in deter- 
mining the ED 2 angular response and pointing. The center of this measured response 
function is 29 deg from the instrument center line and the FWHM dimension is 28 deg. 

The calibration measurement for EDI was slightly better and is shown in Figure 
6. This plot contains both data from EDI and ED 2. The FWHM dimension of the EDI 
curve is 38 deg and the center line of the response is at 11 deg from the instrument 
center line. The data of Figure 6 were taken for Head B of the Plasma Composition 
Experiment and are essentially the sum-total of the preflight, low-energy calibration 
data. 


The foregoing first order effort to model the forward part of the ISEE instru- 
ment was an attempt to establish an analytic base line from which to judge the per- 
formance of the instrument vis-a-vis its design parameters and then extrapolate to a 
response for the very low energy Maxwellian plasma of the plasmasphere . While it 
was intellectually interesting, it probably adds only little to our understanding of the 
instrument. If the problem in the comparison with the calibration data is geometrical, 
that is, if ED 2 is shadowed more than intended by internal surfaces, then the 
geometry of the model could be adjusted until the calibration data were reproduced 
and a valid extrapolation could be accomplished. If, on the other hand, the problem 
is one of stray fields in the energy analyzer, then accurate modeling and extrapola- 
tion is unrealistic. (A. Ghielmetti suggested there should be a deflection in the z 
direction due to fringe fields in the energy analyzer, and E. Shelley adds there may 
be a z component to the field in the edge regions of the RPA-Accelerator section.) 
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Figure 6. Detector angular response calibration, Head B, 
energy detectors 1 and 2. 


Note that the calibration data were taken at 19 electron volts. With lower 
energy plasma the effects of strong fields would be magnified in an unknown fashion. 
Finally, there is the possibility that the calibration itself was in error. It was done 
in the chamber at Bern which, at the time, was apparently experiencing difficulties 
with low energy ions. There was no independent beam monitor, nor was the measure- 
ment checked for repeatability. 

In the case of the second order (ray tracing) analysis, it is difficult to argue 
that the effort should not be pursued. However, it is our opinion, that without 
supporting calibration data the final results of the effort would be subject to the 
same objections as the previously described model. The exercise would be intellec- 
tually interesting but there would be no way of verifying the final conclusions. 

The final recourse is to study the flight data. Comparison of the alignment of 
the peaks of (presumed) trapped distributions with the 90 deg pitch angle crossings 
in the Ram Angle Curves suggests a displacement of about 10 deg between MD and 
ED2 (J. L. Horwitz, private communications). 

Examination of the Ram Angle Curves in regions of the plasmasphere where the 
instrument is switching between EDI and ED 2, shows evidence of a separation in 
pointing between the two detectors. Since the shape of the curves is proportional to 
M sin 0s, it may be possible to derive a reliable &0 between the two detectors; 
although sheath effects could influence this and the picture is complicated somewhat 
if the two detectors do not have the same 0 angle width. 
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For the present we assume that EDI and ED 2 are colinear and are displaced 
from MD by 10 deg. When the sheath influence on the Ram Angle Curves is better 
understood [1], and routine curve fits are possible [2], it should be possible to 
look at the switching data in the plasmasphere as a final verification. (Hugh 
Comfort's analysis of the ISEE data is summarized in his thin sheath paper [2]. The 
influences of external effects [i.e., sheath /charging] apparently were never fully 
determined . ) 


II. INTER-SENSOR GAINS 


Figure 7 shows the approach used for comparing the gain between the various 
detectors on each head. If the spacecraft is in a predominantly H + plasma, and if it 
is switching gain, and if the mass analyzer is set to H + , then it would be possible 
to observe a consistent set of factors between the detectors which would allow all the 
sensors to be related to one another. In going from ED LO to MD HI, the same total 
gain should be observed no matter which way we go around the loop. The numbers 
shown in the figure are approximately the numbers derived from the preflight measure- 
ments and calibrations. 

NOTE: On the RAM angle plots all measurements have been corrected 

according to the particular detector on which they were taken 
and are all adjusted to the most sensitive detector, MD-H1. 

To verify these preflight numbers on flight data, regions in the data were 
located where the instrument was changing ranges (generally near the plasmapause) 
and the plot program was run with gain factors of unity input. Figures 8a and 8b are 
examples of the results. (The RAM angle is off in Figure 8a.) 



•DERIVED FROM THE OTHER THREE 


Figure 7. Summary of relative detector efficiencies for Head B, 
based on pre-flight calibration data. 
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ISEE PCE HEAD B (HYDROGEN) 
DAY 331 OF 1977 
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RAM ANGLE (DEGREES) 


ire 8a. 


Flight data from Head B , comparing data from energy 
detectors 1 and 2. 


ISEE PCE HEAD (HELIUM) 

DAY 331 OF 19/7 
1036 UT R = 2.9 LT = 1306 
L = 3.1 



RAM ANGLE (DEGREES) 


10 


Figure 8b. Flight data from Head B, comparing data from the mass detector 

at low gain and high gain. 





The data taken at the two different settings at zero RPA potential can be 
clearly seen in the plots. The plot labeled H + is the ED detector and the plot 
labeled HE + is the MD detector set on helium. In the two cases a gain of 31.4 is 
measured between ED-HI and ED-LO and 28.2 between MD-HI and MD-LO. A number 
of passes were done and the results are listed in Table 1. (NOTE: There was not 
a great deal of care taken in removing background counts from the data so some of 
the particularly low measurements may be influenced by it.) (This can be very 
important in the plasmasphere . ) 

It is apparent that there is a substantial difficulty in the ED-HI and ED-LO 
comparison. In view of the results of the previous section, the most likely possi- 
bility is that there is a drop in sensitivity in ED-HI. (This was verified at LMSC 
and Lindau.) 

To cross check to the mass analyzer, MD ; data from a special calibration called 
a "bias sequence" was obtained from Day 314, 1977 - outbound. This was a difficult 
section to analyze because it was taken near the plasmapause and in portions of it 
there was a fairly high background count. However, it did support a considerably 
higher ratio between ED-HI and MD-HI, a ratio of the order of 16. Hence, for the 
first pass through Head B data the gains of Figure 7 were replaced by those shown 
in Figure 9. (It is noted that they are not self-consistent but neither were the data) 

In the latter part of the HEAD B data (the last six weeks), a cold plasma mode 
was used in which a few spins of H + were obtained at about 30-min intervals in MD . 

A portion of this data has been studied, three samples of which are shown in Figures 
10a, 10b, and 10c. These are plots of counts versus RAM angle. The top panel 
shows the calculated gain factor and markers at the 90- deg pitch angle locations. 

A factor of 2 seems to be coming out of this set; a number more in agreement with 
the preflight data, although it is likely that a look at a large amount of data will 
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Figure 9. Summary of relative detector efficiencies for Head B, as modified to 
include results from analysis of flight data, early in the detector life. 
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ISEE PCE HEAD 8 
DAY 100 OF 1978 
0951 UT R = 5.7 LT = 2017 
L = 7.2 



Figure 10a. Pitch angle distributions for H+, showing the ratio of the fluxes 
for the mass detector and energy detector 2. 

ISEE PCE HEAD B - HYDROGEN 
DAY 100 OF 1978 
1016 UT R « 6.5 LT = 2037 
L = 8.1 



Figure 10b. Pitch angle distributions for Head B, H+, ED2 and MD . 
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I SEE PCE HEAD B - HYDROGEN 
DAY 100 OF 1978 
1039 UT R = 7.2 LT = 2052 
L = 8.8 



Figure 10c. Pitch angle distributions for Head B, H+, ED2 and MD . 

produce results as variable as did the checks between ED-HI and ED-LO. If the 
results of Figure 10 are to be believed, as well as the results of Table 1, then the 
gain pattern of Figure 11 would appear to be the proper one. 

This, however, brings up one of two very unlikely possibilities: 

1) Detectors MD-HI, MD-LO, and ED-HI have all dropped in sensitivity by 
something close to the same amount . 

2) ED-LO has increased in sensitivity by a factor of 4 or so. 

There is another possibility. It should be noted that the data of Figure 10 
were taken in a region of field aligned flows where the energy of the ambient par- 
ticles was substantially higher than that of the plasma near the plasmapause of Table 
1. Further, at times in the magnetosheath the instrument is observed to auto-range 
between ED-HI and ED-LO. At these times the gain factor appears to be more like 
that of Figure 7. It is, therefore, quite possible that the gain factor between ED-HI 
and ED-LO is energy sensitive by a factor of 6 over the energy range of 0 to 100 V. 
(E. Shelley notes that this is quite possible. Also, he expects the MD-LO sensitivity 
to change (drop) with time more than the others. The change will be mass and 
energy dependent since it is related to the secondary electron production ratio of the. 
incident ions.) 


Tentative conclusions regarding ED-2, Head B, are: 
1) Its acceptance angle is unknown. 
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TABLE 1. GAIN COMPARISONS HEAD B 

DAY OF ED MD 

YEAR GAIN GAIN 

RATIO RATIO 


1977 


1978 


314 
321 
331 
333 
341 
343 
345 
352 
357 
355 
360 
364 
002 
009 
016 
019 
021 
028 
031 
033 
040 
04 2 
054 
059 
062 
064 
066 
086 
100 

MEAN = 
STANDARD DEVIATION = 



16.4 

22.7 

43.3 

31.4 

31.4 

28.2 

43.3 

31.4 

48.3 

31.4 

48.3 

38.9 

48.3 

16.4 

43.0 

43.0 

41.0 

18.3 

43.0 

43.0 

41.0 

22.7 

48.3 

34.9 

34.9 

— 

14.8 

31.4 

10.7 

31.4 

66.7 

10.7 . 

78.4 

13.3 

48.3 

31.5 

10 . 1 

10.1 

29.7 

34.9 

74.2 

15.6 

48.3 

17.4 

74.3 

43.3 

53.8 

7.0 

10.7 

31.4 

14.8 

22 . 7 

12.6 

6.6 

— 

9.6 

25.3 

18.3 

39.4 

24.9 

20.0 

11.3 





Figure 11. Summary of relative detector efficiencies for Head B 
based on flight data, late in detector life. 







2) Its pointing direction is uncertain. 

3) Its relative sensitivity is variable. 

III. DENSITY CALCULATIONS 


The counts/second from a limited aperture RPA set to zero volts retarding 
potential is given by (planar approximation) 


counts/second = 


A t N f 

U) 3/2 a 3 J 


” + * D + 1 71 ‘ 6 o + 1 


/ 


g(e,4>) 


t r + 


D 2 


7i - e - 

o 2 


f 3 ( 

I v exp I - 


a - -[V 2 + V 2 + 2 V Vf ( e,4>)]’ 

x | v exp [ 2 1 dVdedcj) 

a z 


') 


if the ambient medium is a Maxwellian distribution, and 


A = areas of the detector 

t = transmission of the grids 

N = ambient density 

V = spacecraft velocity 
s 

a = mean thermal velocity 

2 2 

g( 0 , d> ) = sin0 cos* cosebsin 8 + sin© sind> sind>sin 0 + cos0 cos0sin0 
b o o o o o 

f ( 0 . <t> ) = sin0 sin0coscj) + cos0 cos0 
s s 

and the several angles are defined in Figure 3. 

This integral has the pleasant property that its value is relatively constant over 
a wide range of Mach numbers (M = V la) . 

Therefore, if the above equation is written 


counts/second = — \ — V {INT} 

L s 
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where considerable rearranging- has been done, it is possible to evaluate the .(quantity 
{INTI at some appropriate point and then use this expression over a considerable 
range at temperatures and flow velocities with little error. The point cKosen gave a 
value of 


UNT} = 0.50 


• • ' • m. »s. I • i • {> : - 1 < , • 

and the error resulting from this approximation is shown in Figure 12 for a range of 
V, T parameter space. 

Another pleasant occurrence is that within the plasmasp here the only perturba- 
tion to the flow velocity of the plasma is from co-rotation flow , or. something close to 
it, and this tends to be smail with respect to the sp.ac.ecraft velocity and also tends 
to be at right angles to it. Hence, the co-rotation velocity tends to shift the incom- 
ing angle of the flow but contribute little to the resultant speed. Therefore, it is 
also a fairly good approximation within the plasmasphere to assume the plasma to be 
flowing at the spacecraft velocity. 

Given these two approximations, which together contribute probably less than 
20 percent error, it is possible to explore the absolute sensitivity of the ISEE instru- 
ment with the expression 

„ _ 4 (counts/sec) 

(At) V c 

b 



Figure 12. Effects of velocity and temperature on the density calculation. 


.16 




If N can be determined from an independent measurement, the quantity, A^. - 

(At) can be derived from the data and compared with preflight measurements. This 
was done for time points in the early data using the upper hybrid resonance measure- 
ments and the results are listed in Table 2. All measurements are for EDI, Head B. 


TABLE 2. EFFECTIVE AREA CALIBRATION 


1977 

TIME 

COUNTS 

DENSITY 

VS. 

A eff 

DOY 

(UT) 

PER SEC. 

(IONS/CM 3 ) 

(M/SEC) 

(mm 2 ) 

329 

? 

29665 

1985 

6931 

8.62 x 10“ 3 

333 

19:53 

14710 

1117 

6015 

8.76 x 10 -3 

336 

5:28 

14719 

1117 

6861 

7.68 x 10~ 3 

350 

13:36 

17994 

1117 

6861 

9.39 x 10 -3 

362 

12:25 

16291 

1117 

7013 

8.32 x 10~ 3 


A » 8.55 x 10" 3 

ef f 


The preflight information was as follows: 

A eff = <P h y sical area °f detector) x (grid transmission) ( number grids) 
x (channeltron efficiency) 

= (30 x 10" 3 mm 2 ) (0.9) 3 (0.33) 

= 7.22 x 10 3 mm 2 


A close comparison! 

Alternatively, this approach can be used to estimate density directly from the 
RAM angle plots if the plasma meets the criteria of corotating, thermal, and Max- 
wellian. The numbers on the plots are actually counts per frame but have been 
adjusted for the relative gain between the detectors such that they appear to be 
coming from MD2. This leads to the equation for density of 


/ counts) 2 

\ At ) (A eff ) (G x ) (G 2 ) {INT} V g 
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where 


At is the time interval of the frame accumulation (0.219 sec) 

is the gain factor between ED 2 and MD2 (G^ = 16 in most of the plots) 

G g is the gain factor between EDI and ED 2 (G2 = 47 for most of the plots). 

{INT} = 0.50 


N = (3.36) counts 


3 

where V is in meters/sec, N is in ions/cm , and counts are from the plots. 

O 

This approach was used to determine the density in the plasmasphere for 
several passes and the results were compared with the wave data by Hugh Comfort. 
An example of the results is shown in Figure 13. 
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Figure 13. Comparison of plasmasphere density profiles obtained from 
the PCE and the plasma wave instrument. 



This particular pass is admittedly one of the better comparisons, however, they 
all at least were "in the ball park." In general, the. results showed that the Plasma 
Composition Experiment was overestimating the density by about a factor of two when 
ambient densities were in the 1500 to 2000 ions /cm 3 range and underestimating the 
density by a factor of perhaps 4 or 5 when the ambient density was in the 100 ion /cm 3 
range. The error is at least in the right direction to be explained by spacecraft 
potential. At present the best guess is that the spacecraft is of the order of 0.5 V 
negative at the higher density range at 1000 to 2000 ions/cm 3, increases to the 
neutral point at about 100 to 50 ions /cm 3 and goes positive at smaller densities. This 
question will presumably be answered if Mozer's data is ever figured out. 

Singh has looked at the enhancement of current to a limited aperture RPA when 
the spacecraft is negative. Table 3 lists some preliminary results for a 45 x 45 deg 
aperture. They are given for comparison purposes. 


TABLE 3. SATELLITE POTENTIAL EFFECTS ON CURRENT TO AN RPA 


q<J> s /kT 

-1 

-2 

-3 

-4 

-5 

-10 

Enhancement 

2 

3 

3.9 

4.7 

5.5 

8.3 


Of course, from physical considerations it would be expected that the crossover 
in density determinations would occur at the point where the spacecraft goes through 
zero potential. If there were to be forced on the data, it would tend to lower the 
effective area of EDI by about a factor of four and would indicate enhancements in 
the density measurements by as much as a factor of 10 in the vicinity of 2000 ions/ 
cm 3 . 


Finally, Table 4 is from a note obtained from W. K. Peterson. These numbers 
give his last guess for geometric factors for MD and ED for Head A and Head B for 
several mass species in the low energy channel of the instrument. These were 
apparently derived from a combination of calibration fits to analytic instrument models 
and were influenced by flight data. They are for the low gain detectors. The 
analysis assumes the energy passband is filled. The ratio of GSE for MD and ED 
should be observed in the data if the environment is all H + and MD is set on 57. 
These data are the result of best fits for 5 or 6 magnetosheath intervals in October 
and November of 1977, counts/second = G<5E flux/keV. 

TABLE 4. GEOMETRIC FACTORS 




For ESTEP = 0 




HEAD A 

HEAD B 

ED 


Ec = 49 eV 

Ec = 54 eV 


GSE 

0.19 x 10' 3 

0.20 x 10-3 

MD 


Ec = 40 eV 

Ec = 44 eV 

H + 

GSE 

0.36 x 10'i? 

0.37 x 10“ J 

He ++ 

GSE 

0.45 x 10, 

0.53 x 10, 

He + 

GSE 

0.72 x 10, 

0.56 x 10, 

0 + 

GSE 

0.31 x 10' d 

0.52 x 10' d 
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IV. ENERGY PASSBAND-BIAS CHECKS 


There is an influence of the instrument on the energy distribution of incoming 
particles in the lowest channel of the electrostatic analyzer. This channel, nominally 
100 eV in width but more likely about 150 eV in width, has a response versus energy 
which is Gaussian. Although the influences of this characteristic are probably only 
important for high temperature, low Mach number plasma, it should be noted. The 
only measurements of this response were done in flight during the "bias sequences" 
mentioned previously. These measurements were done within the plasmaspher'e in. the 
presence of a cold plasma and the bias on the electrostatic analyzer was moved in 
25- V steps to sweep the passband across the incoming cold beam. The measurements 
from DOY 314, 1977, have been analyzed. 

Figure 14 is an example of the results of one such sequence (J . L. Horwitz, 
private communication). This one is for MD Low Gain, Head B. The x-axis is steps 
(approximately 25 V), the y-axis is counts. 

The other possible area of influence of the energy distribution is from sheath 
effects. These are currently under study by Singh and would be expected to be 
small at the higher energies but could be important at the lower energies of the 
plasmasphere, particularly in the area of determining the spacecraft potential from 
the ion data. 


10.000 
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4.000 

3.000 

2.000 
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2866 2891 2915 2934 2965 2990 3015 3033 BIAS VOLTAGE 

Figure 14. Characteristics of the electrostatic analyzer response in the 

RPA mode, Head B, flight data. 


ENERGY PASSBAND OF THE '0-100 eV' CHANNEL 


MASS DETECTOR - LO GAIN 
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V. HEAD A CALIBRATION 


(Predicted Energy Response of the Sharp ISEE Instrument 
at the Lowest Energy Step) 


Figures 15 and 16 show the predicted ED2 and MD energy response at E-STEP 
00000 for sensors A and B, respectively. The vertical scale is counts normalized to 1; 
the horizontal scale, energy in eV, is in units of 20 eV per major division. Rather 
than fix the external energy scale and draw eight response curves, one for each 
B-STEP, one response curve requiring that the energy scale be shifted appropriately 
for each B-STEP has been drawn. The eight points indicated on the abscissa are the 
zero points corresponding to the eight subscripted B- STEPs. Thus, for sensor A, 
the ED2 energy response centerline is +13 eV at B-STEP 110, -4 eV at B-STEP 111, 
etc. The ED 2 and mass analyzer (MA) energy response curve shapes are scaled from 
measured data taken with an 2^ + beam at E-STEP 00001 and B-STEP 011 (sensor B) 
or B-STEP 101 or 100 (sensor A). 

The scaling was based on the following equations: 


E^ n = E^ + V c = constant at any given E-STEP 


( 1 ) 


E 0X E 1X 


\ AV l ) c V AVj / 


( 2 ) 


V c = V(R C ) = AV 


ln ( R c /R inner ) 
ln ^ R outer //R inner^ . 


- V. - AB 

inner 


(3) 


K. = 


E. . E . + V . 

in c l xci ci 


i AV. 

i 


AV. 

l 


(4) 


where 

E = particle energy 

V = voltage or energy depending on context 
A V = potential difference between ESA deflection plates 
AB = voltage or energy difference between B-STEP settings 
R = radius of ESA plates or center trajectory 
K = analyzer constant 
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21-OCT— 77 EH 
2H + 23— JUNE— 77 



B ENERGY (eV) 

Figure 15. Characteristics of the electrostatic analyzer response in the 
RPA mode, Head A, ED2 and MD, laboratory. 



Figure 16. Characteristics of the electrostatic analyzer response in the 
RPA mode, Head B, ED2 and MD, laboratory. 





and subscripts 

in = internal 
X = external 
c = center 
i = E-STEP index 
0 = E-STEP 00000 


1 = E-STEP 00001 
inner = inner ESA plate 
outer = outer ESA plate. 

Particles entering on a centerline trajectory are assumed to be preaccelei’ated 
to the center potential in the ESA, V . Inner and outer ESA plate voltages and 

AB's were measured in atmosphere at 10°C prior to calibration in vacuum. These 

were used to calculate AV. and V . E 1V was taken to be the center of the full 

i c IX 

width at half maximum (FWHM) of the EA energy data measured in vacuum. Energy 
values for those data were read from a digital voltmeter indicating source potential 

of the beam. Initially the analyzer constant K of each EA and M A was assumed 

constant for all i. This allowed E QX to be calculated from E^ x at the B-STEP at 

which E 1X was measured by applying equation (2). FWHM Q = (^oin^lin^ FWHM^, 

approximately 6.5 percent narrower. Finally, the zero positions on the scale for all 
B- STEPS are determined by adding or subtracting the appropriate AB. 


It was discovered during calibration in vacuum with a particle beam, however, 
that the high-voltage power supplies were not quite sufficiently temperature com- 
pensated, and thus during operational warmup AV. diminishes. This problem was 

treated using equation (4) as though AV^ were constant and variable. The percent 

variation in K. is a measure of percent variation in 4V.. versus E and E-STEP is 

plotted for sensor A, EA and MA in Figures 17 and 18, respectively, and for sensor 
B in Figures 19 and 20, respectively. The range of variation at the E-STEP 00000 
position was assumed equal to the 00001 data range plus 25 percent margin on either 
side for sensor A and 50 percent margin on either side for sensor B . This range of 
uncertainty is reflected in Figures 15 and 16 by the sketches in the upper left-hand 
corner depicting a band of uncertainty around any given zero point. The arrow 
indicates the direction in which the zero tends to shift as the instrument warms up . 
For sensor A, the uncertainty band is explained thus. 


-DIRECTION 0 POINT TENDS TO SHIFT AS 
INSTRUMENT WARMS UP 

— 25% SAFETY MARGIN, BOTH SIDES 

— FULL RANGE OF Kl DATA 


For the initial turn-on, B -STEPs 101 for sensor A and 011 for sensor B are proposed. 
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Figure 19. Analyzer constant for Head B, energy detector (ED) 
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Figure 20. Analyzer constant for Head B, mass detector (MD) 





For sensor B, the EDI energy response was determined by fixing the beam 

energy at 19 ± 1 eV and measuring count rate as a function of B-STEP. The experi- 
mental EDI data and the predicted ED 2 response curve are plotted in Figure 21. 

They suggest that the EDI response peak is about 21V higher than that of the ED 2. 
There are no other calibration data against which to check this result. The EDI 
response is represented on Figure 16 by scaling the ED 2 response curve to the 
correct width and location. 



Figure 21. Characteristics of the electrostatic analyzer response in 
the RPA mode, Head B, EDI, laboratory. 


VI. RAY TRACING 


The ray tracing analysis for the PCE is as follows: 

1) Refer to Figures 22, 23, and 24 for diagrams and dimensions. 

2) 0^ is measured in the x-z plane and e 2 is measured in y-z plane. 

3) Assume that E = E , E =E = 0 between B2 and G, between G and B3, 

z x y 

and between B3 and SI. Assume E = 0 between the entrance aperture and B2, 
between SI and the analyzer entrance, and between the analyzer exit and S2. 

4) Assume that the Fringing Field correction at the analyzer entrance and exit 
can be made by a refraction approximation. 
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Detailed Procedure 

1) Let x^, y ^ , be coordinates in the plane of the entrance aperture and let 

x^, y 2 k e coordinates in the plane of grid B2. To select a particle, choose at 
random x^, y^, x^, y 2 subject to the restrictions -4.72 = x^ =4.72; -1.1 - =1.1; 

-3.7 - x 2 ^ 3.05; and -0.263 ^ y 2 ^ 0.263. 

2) Let 0 y be the angle of the particle direction in the x-z plane and 0 2 be the 
angle of the particle direction in the y-z plane. From the geometry: 



-40.9° = = 38.7° 

- 8 ° - 0 2 = 8 ° . 
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3) Set up 3 matrices with coordinates 0^, 02 in 1° x 1° increments. M q will 
record the "unsuccessful" particles, will record the "successful" particles that 
strike EDI and will record the "successful" particles that strike ED 2. 

4) To start the computation, find the initial location and velocity components 
of grid B2. These are: 


z o = - 3 - 0 


x = x„ 
o 2 


y o = y 2 ‘ 


Define D = [9.7 2 + (x x - x 2 ) 2 + (y.^ - y 2 ) 2 ] 2 . 
between the entrance aperture and B2). 


(D is the total distance traveled 


V q = initial particle velocity 




V 


xo 


V o (X 2 • x l } 
D 


-V sin0 1 
o 1 


V o <y 2 ' y l> 


y° 


D = ' V o Sin0 2 


CO 


■ v o (¥) 


5) 


Between B2 and G, 



-V 

E = - RPA 
z Az x 299.7295 


-V 

RPA 

149.86475 


The equations of motion are: 


a) 

x' 

= X + V 1 
0 XO 

b) 

y' 

= y ^ 

+ V 



0 

yo 

c) 

z' 

= Z + V t 

o zo 
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d) 


vxo 


e) v , = vyo 

f) V , = vzo + 

z' 

z' = z(G) = -2.5; Z Q = z(B2) = -3.0 . 

Equation C can be solved for t giving; 


-V 


r ± /V ' 
zo zo 


+ 2 ( Z 1 


Z ) — E 
o m z 


eE 


m 


2 

Case A : b - 4 ac < 0 -* no real solution for t exists, i.e., the particle is 

stopped by the retarding potential . 

2 2 1 

Case B: b - 4 ac > 0 and [b - 4 ac] 2 - | v | there will be two positive 

zo 

solutions for t. Pick the smallest positive value. 

2 2 1 i i 

Case C : b - 4 ac > and b - 4 ac 2 > |v zo | -+ there will be one positive and 

one negative solution for t. Pick the positive solution. 

After finding t (assuming a real value exists), use equations (a), (b), and 

(f) to compute x', y', v '. 

z 

If v Rp = 0, then the quadratic solution is undefined. In that case use: 


zo 


and then use equations (a), (b), and (f) to compute x\ y\ v '. 

z 

6) At grid G, a new computation is started. So redefine the initial positions 
and velocities according to: 


x = 
o 


x' 


= y’ 


V 

xo 


V 


X 



From results of part 5 
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From results of part 5 


V yo v y' 


V = V , 

zo z' 


z' = -2.0 


z = -2.5 
o 


E 

x 



E 

z 


V 

RPA 

Az x 299.7295 


V 

RPA 

149.86475 


Repeat the procedure outlined in (5) above to trace the particle from G to B3. 

Again, pick the smallest positive real value for t. B3 is a limiting aperture, and so 
at B3 need to check if -3.7 < x' < 2.88 and -0.1625 < y < 0.1625. If both these 
conditions are not met, then the particle was not "successful." 

7) At Grid B3, a new computation is started. So redefine the initial positions 
and velocities according to: 


x = 
o 


y 


o 


= y' 


V 


xo 



V 

yo 



V = V 1 

zo z 


From results of part 6 


z’ = -1.0 


2935 E.G.S. 6/23/76 


z = -2.0 
o 


E = E = 0; E 
x y z 


3000 


299.7295 


= 10.0069 statvolts/cm 


Trace the particle from B3 to SI using the procedure outlined in (5) above. Since 
this is an accelerating potential region, the quadratic solution for t will result in one 
positive and one negative solution for t. Pick the positive solution. SI is a limiting 
aperture, and so at SI need to check if -3.7 < x f < 2.765 and -0.075 < y' < 0.075. 

If these conditions are not met, the particle is not "successful." 

8) Between SI and the analyzer plate entrance, the particle is in a zero-field 
region. Re- initialize the position and velocity according to: 
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The equations of motion are: 


a) 

x' = 

X + V t 

o xo 

b) 

y' = 

0 

>> 

> 

+ 

o 

>> 

c) 

z' = 

Z + V t 

o zo 

d) 

V ' 

z 

II 

< 

tS) 

o 


From results of part 7 


Solve (c) for t and then use (a) and (b) to find x', y'. Also, v ' = v , v ' = v 

J x xo y yo 

9) Tr acing Particle Through the Analyzer Plates 

At the entrance to the plates, change from cartesian to cylindrical coordinates 
(r, <J>, O . The electric field is radial, = -A/r, E<j> = = 0. 

1 / V 2 V 1 \ _ 2Vdef 
r r \ln(r 2 /r 1 )/ r lntr^r^ 


The potential is: 
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The potential outside of the plates is <p Q - -3kV = -10.0069 statvolts. The 
potential inside the plates is <f>(r) which is not in general equal to <j> . Hence, there 

is a refraction at both the entrance and exit of the plates which causes a change in 

v and hence a change in angle. Let x,y,z=0,v , v , v be the initial 

zo & 5 o^oo xo yo zo 

coordinates and velocities in the analyzer plate entrance. The effect of refraction is 

to change v ^ to v zo ' according to the equation: 

V ' = I"-- (<j> - <Kr)) + V 2 1 2 

zo |m T o T zo J 

where 


r = 7. 27 + y 
°o 

<p = -10.0069 
o 


4> (r) 


-11.42 In hH- 8.6007 In JL* 
r b.77 

0.137769077 


Now transform the initial conditions to cylindrical coordinates according to the 
equations: 


a) r Q = initial particle radius = 7.27 + y Q 


,2 A „ 2 ^ 


b) v q = initial particles velocity in the r - <j> plane = (v zQ ' + v^ o ) 


c) v = v 

Co xo 


d) c = x 
s o o 

-1 V o 

e) a = initial angle of the v vector relative to the z axis, a = tan 

n & o o V 

zo 


The following is based on the treatment of Hughes and Rojansky [3], The 
differential equation describing the orbit is: 


. 2 2 /d 2 u \ 
h a I — 7) + a I 

\ d<j> / 


Ae 

m 


y= F 


h = r v cos a = initial angular momentum 
o o o & 
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Making the substitution y = m/p 


= =_! . y 

d, 2 » y 


o 2 = 


Ae 


2 2 
mv cos a 
o o 


with the boundary conditions: 


<+3 = 0, y 



tan a 

o 


The equation has the approximate solution: 

g) y = c + (1-e) cos /2 <p - (tan a Q sin /2 <}>) / /2 

h) ^ .= - tan a = -(1-c) /2 sin /2 <j> - tan a cos /2 <j> 

a 9 o 

At the exit of the plates, <f> = 111.7° 

cos /2 <j> = -0.926972 
sin S2 <p = 0.375130 

Equations (g) and (h) are used to compute y^ and hence r^ = r 0 ^f anc * a f = '^ an * 
(dy/d<{))n at the analyzer exit. Since we assumed that E 5 0, then V „ = V 

Now we have Vp r^, a^, V £. Define a new cartesian coordinate system at the 
analyzer exit with Z = 0 at the exit plane. The transformation is: 


V = V cos a f 
zo o t 


V = V sin a,, 
yo o f 
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V 

xo 



X = X (at entrance) + V t 
o o p 


Y = r - 7.27 
o 

Z = 0 
o 


The transit time t 

P 


through the plates can be approximated by: 


t 


P 


path length 
V 

o 


14.173 

V 

o 


As the particle leaves the plates, there is a refraction due to the jump in potential 
from 4>(r^) to <j>j = -10.0069 statvolts. Hence as before, compute a new V zQ according 
to: 

V ’ = f — (<Kr f ) - <j> n ) + V 2 1 2 
zo I m r ^1' zo J 


10) Between the exit of the plates and the plane of EDI and ED 2, the particle 
is in a field- free region. The equations of motion are: 


a) 


= X 

+ V 

F 

o 

xo 

b) 

Yp 

= Y 

+ V 

F 

o 

yo 

c) 

Zp 

= Z 

+ V 

F 

o 

zo 

d) 

Z 

o 

= o; Z F = i 


Solve equation (c) for t and then use (a) and (b) to compute Xp and Yp, the final 
position of the particle in the plane of EDI and ED 2. 

11) For EDI, the condition for a "successful" particle is that it strikes within 
a rectangular aperture of Ax = 0.2 cm, Ay = 0.5 cm centered at x = -2.79, y = 0. 
For ED2, the condition for a "successful" particle is that it strikes within a 0.319 cm 
diameter circle centered at x = -3.7 cm, y = o. 
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APPENDIX A 


DEFINITIONS AND PHYSICAL CONSTANTS 

E = electric field, c.g.s. units, statvolts/cm 
V = particle velocity (cm/sec) 

= initial particle energy (ergs) 

V q = initial particle velocity = 

(To convert from energy in electron volts to energy in ergs multiply by 1.6022 x 

nr 12 .) 

- 24 

M = proton mass = 1.673 x 10 grams 

e = electronic charge = 4.80325 x 10 ^ e.s.u. 

1 statvolt = 299.7925 volts 

Section 9 VdeF = 423.6399 volts 

For E Center = 3075 eV 3000 E.G.S. 6/23/76 



A 



20.51419 


2 _ Ae _ 5.889706 x 10 15 

C „ 2 2 „ 2 2 
m V cos a V cos a 

o o o o 


preceding page blank not note 

V/ 


^ GS ±2JL^mmiONALL)g BUNK 
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APPENDIX B 


PARTICLE "SCORECARD” MATRICES 


-40.6 ^ 6 1 £ 3.87 
-8° i e 2 ^ 8° 

1. Compute 0^, 0£ . 

2. Round off and fix according to: 


e 

e 


1 

2 


= FIX( 0^ 
= FIX(0 2 


+ 41°) 

0 ^ 

e 1 $ 79 

+ 8°) 

0 ^ 

0 2 $ 15 


Define M.. 

i] 

i = 0 i + 1 1 $ i ^ 80 

y=0 2 +l 1 £ ] £ 16 

Then the 0 windows are computed from the i, j values according to: 


1L 

= (i - 

1) 

- 41 = i - 42 

1U 

= (i - 

1) 

- 40 = i - 41 

2L 

= (j - 

1) 

a> 

i 

ll 

00 

2U 

= (j - 

1) 

- 7 = j 1 - 8 


Determine the probability distribution of 0^ 2 in ISEE ray tracing. 
. (x. - x„) 


0^ = tan 


9.7 


Xj is a random variable uniformly distributed over the interval - 472 £ X^ 
X <2 is a random variable uniformly distributed over the interval - 3.05.$ X 2 


4.72. 

< 3.05. 


PRECEDING PAGE BLANK NOT FB MS) 
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Define y = - Xg. 

What is p(y)? 

p(y) can be found by a graphical technique discussed by Pfeiffer [4], 
Chapter 3. 



The probability density function p(y) is proportional to the length of the line enclosed 
within the probability mass rectangle. 

X 1 = Y - 3.05 

X 2 = -Y + 4.72 

For 1.67 ^ y ^ 7.77 

L = [ (x x - 4.72) 2 + ( x g - (-3.05))V 

L = /2 (7.77 - y) 1.67 < y < 7.77 

L = (6.10) 0 < y < 1.67 
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For 


p( e i ) = = P(y) 9.7 (1 + tan\) 

p(9) = P(y = 9.7 tan 0^ 9.7 (1 + tan 2 0 1 > 
3.87° < 8 < 38.7° 
pCep = 0 for 1 8jJ > 38.7° 


In a similar manner we can compute p(0 9 ). 





Define Z = Y 1 - Y 2 . 

Y 1 is r.v. uniform over -1.1 ^ y ^ 1.1 
Y 2 is r.v. uniform over -2.63 ^ y 2 i 0.263 



for 0.837 ^ z ^ 1.363 

(-0.263)) 2 ]* 

1.363)]* 

For 0.837 < z < 1.363 
For 0 < z £ 0.037 

0.837 £ z £ 1.363 
0 ^ z ^ 0.837 


y x = z - 0.263 


y 2 = -z + 1.1 


L = [( yi - l.ir + (y ( 


[ ( z - 1.363) z + ( - z + 


L = /2 ( 1.363 - z) 


L = /2 (0.526) 


, N L 1 
P(z) = ^ ri572 


p(z) = 0.864 (1.363 - z) 
= 0.455 
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APPENDIX C 


COMPUTATION OF GEOMETRIC FACTOR (G.F.) FROM THE DATA 


G. F 


d l F y*A(fi)dfi = Ja(Q 


1’ 


e 2 ) 


do., d0 f 

1 C 


In the ray- tracing method used here, the detector aperture is effectively the 
grid B2. The entrance aperture serves only to limit the incoming particles in angle. 

If we approximate the particles in a 1° x 1° block in 0^, e 2 space as a uni- 
directional beam, then 


A(0 r 


0 2 ) 


counting rate (#/ sec) 
flux (#/cm 3 - sec) 


Define: 

= Number of successful particles at angle 6^., 0^. . 

T.. = Number of trials at angle 0.., 0 9 - . 

1] lx 

Now since B2 is uniformly illuminated with particles at constant 0^, 0 2 (X 2> Y 2 has 
uniform joint probability mass over B2), then 


T.. 

"F " ( #/cm^l = 

ij ( ' ff/ ; A(B 2) 


"counting rate"(#) = 


Therefore, 

G.F. 


Z H . . 

— i 


rr t.. 

i.i i) 


G,F * 57.3 X 57.3 X 3-49 5 T.! 

GJ i) 


= 1.0631 x 10 


o H •• 
-3 > _i] 

~ T.. 
i>] ij 
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APPENDIX D 


TRAJECTORY THROUGH PLATES 


Cylindrical Electric Field 



E 


r 




By Gauss' Thm. 


E 


r 


K 

r 



dr = 


v 


2 



v 


2 


The DE of 


h 


2 2 
M 


the orbit is: 



M 


where y 


1/r; h = r Q v o cosa, A 



This can be transformed into: 



Y = 



where C 


Ae 



rav 

o 


where the b.c. at <j> = 0, y = 1, dy/d$ = -tana 

This equation has the approximate solution: 
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INTENTIONALLY BUM 


y 


C + (1 - C) cos ( 2)} 


<p - ( tan a sin /2 <p ) / / 2 


At the entrance, <j> = 0 at y = 1. 
At the exit, = 111.7°, 

cos /24> = -0.92697 
sin m = 0.37513 . 


Now, the condition for a particle to be in a circular orbit is: 



2 

Hence, if C = Ae/MV Q 2, then C < 1 the particle energy exceeds the circular orbit 
energy and C > 1 -*■ the particle energy is less than the circular orbit energy . 

Example: suppose C < 1, a = 0. 

Then y = C - 92697 (1-C) 

U 

= C(1 + 0.92697) - 0.92697 < 1 
But since y = r Q /r, -> r > r Q at the exit. 




So at the analyzer entrance, we know y = 1, r , a. 

At the analyzer exit, we can compute y and hence r and 


^ = -(1-C) /2 sin /2 4> ~ cos /2 <p 

d<p / 2 

which is equal to -tana at the exit. 
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APPENDIX E 


EFFECTS OF REFRACTION AT ANALYZER ENTRANCE AND EXIT 




r 2 r 

v. In — + (v 0 - v.) In — 
1 2 1 r^ 


In 


2 r 

v. In — + v 0 In — 
1 r 2 r, 


In 


v , 


V. 


v 2 


- 3kV 


- 3kV 


V , „ 
def 


- 3kV + V 


def 


V def = 423 ' 64 


-3423.64 volts = - 11.4200 statvolts 
-2576.36 volts = -8.59381 statvolts 

= 10.0069 statvolts 



For particle in region 1 


For particle in region 2 

1/2 mv zl 2 + q^> 1 = 1/2 mv^ 2 + qcj> 2 


-^3kV 
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= 6. 77; r 2 = 7.77 

V. = -3 kV - V . P = -3423.64 volts = - 11.420 statvolts 

1 def 

V„ = -3 kV + V , e = -2576.36 volts = -8.59381 statvolts 

2 def 

-3 kV = -10.0069 statvolts 
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There will be a potential step function at the analyzer entrance and this step function 
affects only the z-component of the velocity, i.e., there is a refraction of the particle 
trajectory. Conservation of energy: 

2 2 
1/2 mv^ + e<j>i = 1/2 mv^ + 

V z 2 = [ 2 m <*1 ' »2> + V zl 2 ]’ 


where 


<j>^ = -10.0069 statvolts 

-11.42 In 1^1- 8.59381 In 
^ 2 = 0.137769077 


Algorithms: 

A. Entrance of Analyzer Plates 

1. Compute x o , y Q , v xq , v zQ at analyzer entrance as before (line 295). 

2. Define r = r Q = 7.27 + y Q . 

3. Compute v according to equation 1 above with v , = v ,v 1 = v 9 

ZO Z X ZO ZO Z Z 

4. Proceed with calculation through the plates as before. 

B . Exit of Analyzer Plates 

1. Compute r^, v zQ , v^ o as before (line 318). 

2. Define r = r^ (line 315). 

3. Compute v * according to equation 1 above with v 1 = v ^ = v 9 . 

ZO Z X ZO Z Z 

4. Proceed as before to trace particles back to detectors. 
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